ADA076597 


OPERATIONS 

RESEARCH 

CENTER 


Du.ri'  •  !;.ur.e;x  n  A"  -  ■ .  ;„:i 

Appwvad  tot  public  ’  ,js  !  t 
Dutnbut’.or  ' 


>0  D  C 

r-srp(7^.  •  ■  iJjd-1-4 

Who*  13  *M\\\ 

n  i  _ t- m— 1 i  \ 


UNIVERSITY  OF  CALIFORNIA  •  BERKELEY 

7  9  11  0  9  0  2  6 


SOLUTIONS  OF  SPHERICAL  INEQUALITIES 


by 

DAVID  GALE 


ORC  79-10 
SEPTEMBER  1979 


i-  i  _ : 

•  .  V 


_ 


SOLUTIONS  OF  SPHERICAL  INEQUALITIES 


by 


David  Gale 

Department  of  Industrial  Engineering 
and  Operations  Research 
University  of  California,  Berkeley 


Acoessicn  Fcr 

NIIS  0.  -..il 
DI'C  TAB 
Unc  .  jounced 
Juci  : : 


Dir.t.-  ’  ' 


SEPTEMBER  1979  ORC  79-10 

This  research  has  been  partially  supported  by  the  National  Science 
Foundation  under  Grant  SOC78-051^6  and  the  Office  of  Naval  Research 
under  Contract  N00014-76-C-013V' with  the  University  of  California. 
Reproduction  in  whole  or  in  part  is  permitted  for  any  purpose  of  the 
United  States  Government. 


<eCu*iTv  classification  OF  This  PAGE  rlFhan  Dm  entered) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

liMTUf.i  IJHI'I  i  j.^— — —— — rriMTlT-f 

mkmmMkJm 

S.  RECIPIENT'S  CATALOG  NUMBER 

&&2±^SSSSSSSSM% 

v  / 

-  — 

7.  AUTHORfAJ  J 

David/Gale  j  S' .  c4  S 

4  -  /  v  — — 

i6  s&WVtfF  I 

s.  performing  ORGANIZATION  name  ano  aooress 

Operations  Research  Center  7 

University  of  California 

Berkeley,  California  94720 

io.  program  element.  project,  task 

ARC  A  *  #0*K  UNIT  NUMBERS 

II.  CONTROLLING  OFFICE  NAME  AND  AOORESS 

National  Science  Foundation  f  It 

1800  G  Street,  N.W.  V. 

Washington,  D.C.  20550  ^ — * 

wmw - Jail 

1*.  MONITORING  AGENCY  NAME  «  AOORESSflf  dllletent  (tom  Controlling  Otllce) 

(t2j\iqbj 

15.  SECURITY  CLASS,  fof  report; 

Unclassified 

1  16.  DISTRIBUTION  ST ATEMEN T  fof  thlt  Report)  1 

Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  ST  ATEMENT  (o(  the  ebetrect  entered  In  Black  20,  II  all  lor  an  I  tram  Report) 


It.  SUPPLEMENTARY  NOTES 

Also  supported  by  the  Office  of  Naval  Research  under  Contract 
N00014-76-C-0134. 

If.  KEY  WORDS  (Continue  an  raaataa  aide  II  neceeeery  and  Identity  by  block  number) 

Quadratic  Programming 
Circumscribed  Ball 


JO.  ABSTRACT  (Continue  on  rarer  re  aide  II  nacaaeory  end  Identity  by  black  number) 

(SEE  ABSTRACT)  ,  ■  ,  f 

I  S  ■  u  'i  ,  1 


7v5~0 


on 

BO  I  JAM  71 


1473  tOITION  OF  I  NOV  SS  IS  OBSOLETI 

t/N  0102-LF -01 44601 


_ Unclassified _ 

SECURITY  CLASSIFICATION  OF  THIS  RACE  (than  Been  Bntered) 


SOLUTIONS  OF  SPHERICAL  INEQUALITIES 


by 

David  Gale 

There  are  many  examples  of  natural  geometrical  problems  which  can  be 
solved  by  linear  or  quadratic  programming.  For  example  the  problem  of 
finding  che  largest  ball  contained  in  a  given  convex  poly tope  (the  inscribed 
ball )  turns  out  to  be  a  linear  program,  while  the  problem  of  finding  the 
smallest  ball  containing  a  convex  polytope  (the  circumscribed  ball )  is 
equivalent  to  a  quadratic  program  [1].  In  this  note  we  extend  the  result  of 
[1]  by  showing  that  the  following  problem  can  be  solved  by  a  simple  quad¬ 
ratic  program: 

Problem: 

Given  a  finite  set  of  balls  in  Rn  decide  whether  their  intersection 
is  nonempty  and  if  so  find  a  point  in  it. 

To  formulate  the  problem  algebraically  let  c^  be  the  center,  the 

radius  of  the  ittl  ball,  i  ■  1,  ...,  n  . 

Problem: 

Solve  the  system  of  inequalities 

2  2^ 

(1)  (x  -  c1)  <  pi  ,  i  -  1,  . . . ,  m  , 

or  show  that  no  solution  exists. 

An  obvious  modification  of  this  problem  is  to  find  a  vector  x  and 
number  z  such  that 

+In  this  paper  the  scalar  product  of  vectors  x  and  y  is  written  xy 

2 

the  scalar  product  of  x  with  itself  is  x  . 


I 


(2) 


z  is  a  maximum 


subject  to 

(3)  (x  -  ct)2  -p2+z<0,i-l,  ...,m. 

Letting  (x,z)  be  the  solution  of  (2),  (3),  it  is  clear  that  (1)  has 
a  solution  if  and  only  if  z  >  0  .  For  the  special  case  p^  ■  0  for  all 
i  it  is  also  clear  that  z  is  negative  (unless  all  c^  are  equal)  in 
which  case,  by  definition,  x  is  the  center  and  -z  is  the  radius  of  the 
circumscribed  ball. 

The  point  of  this  note  is  to  show  that  the  problem  (2) ,  (3)  has  a  dual 
which  is  a  very  simple  quadratic  program  and  hence  can  be  solved  by  well 
known  techniques.  It  then  follows  that  if  the  c^  and  p^  are  rational 
then  (2),  (3)  have  a  rational  solution  which  is  perhaps  not  immediately 
obvious . 

Our  result  is  the  following.  Consider  the  quadratic  program  of  finding 

a  nonnegative  m-vector  X  *  (X, ,  . . . ,  X  )  and  a  number  u  such  that 

1  m 

* 

(2)  u  is  a  minimum 

subject  to 

<3)*  „  -  (i  hsy  ♦  i  -  c\) 

and 


l  Ai  -  1  |[  is  short  for  [  j  . 


it  it 

We  assert  that  (2),  (3)  and  (2)  ,  (3)  are  dual  problems.  More 
precisely: 

Theorem: 

Let  (X,U)  solve  (2)*,  (3)*.  Then  (x,y)  solves  (2),  (3)  where 

*  ■  l  Vi  • 

■We  first  prove  the  usual  inequality  (so  called  "weak  duality").  Let 

* 

(x,z)  be  any  solution  of  (3)  and  (X,u)  any  solution  of  (3)  .  Multiplying 

(3)  by  X^  and  summing  gives 

Oil  Xt(x  -  c4)2  -  l  X^2  +  x  -  x2  -  2  (l  *  l  ^('i  -  »i)  +  * 

’  (*  *  l  Xlcl)2  -  (l  Xlcl)2  '  l  Xi(°l  '  cl)  +  2 

-  (x  -  l  Vi)2  +  2  '  11  i  1  '  U  *  (fr0m  ) 
so 

(4)  Z  <  VJ  . 

_  «.  *  *  —  r* 

Now  suppose  (X,w)  is  a  solution  of  (2)  ,  (3)  ,  and  let  x  ■  l  X^c^  * 
The  Kuhn-Tucker  conditions  for  this  problem  assert  the  existence  of  a  number 
u  such  that 


(5) 


Multiplying  (5)  by 


u  <  2xc. 

2 
c . 

for  all 

—  i 

i 

i 

u  ■  2xc^ 

2 

2 

+  pi- 

Ci 

if  X 

X^  and  summing  gives 


> 


i 

0 


and 


- 


4 


u  -  2i2  +  J  At(p2  -  c2)  , 
or 

(6)  u  -  x2  -  x“  +  l  xi(c^  "  ci)  "  v*  £rom  (3)  • 

-2 

Subtracting  x  from  (5)  gives 

-2  2  -  2 
U  -  X  <  Pj  •  (x  -  Cj) 

so  from  (6) 

(x  -  ci)2  -  p2  +  u  <  0 

so  (x,u)  is  a  solution  of  (3)  and  from  (4)  it  is  optimal,  completing  the 
proof.! 

Remarks : 

*  * 

1.  The  fact  that  problems  (2),  (3)  and  (2)  ,  (3)  are  nonlinear  duals 
can  be  derived  by  appealing  to  general  nonlinear  duality  theorems,  e.g., 

[2]  and  making  some  further  simplifications.  For  this  special  case,  how¬ 
ever,  we  prefer  the  self-contained  proof  which  uses  only  the  standard 
optimality  criterion  (Kuhn-Tucker  conditions)  for  quadratic  programs. 

2.  One  might  hope  to  generalize  the  result  of  this  note  by  replacing 
the  spherical  inequalities  of  (1)  by  general  positive  definite  inequalities 
so  that  the  spheres  of  (1)  become  ellipsoids.  The  following  interesting 
example  due  to  G.  Bergman  shows  that  there  is  no  hope  of  doing  this.  The 
following  inequalities  describe  ellipses  in  the  plane. 
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2  2 

2x  +  3xy  +  2y  -  4x  -  4y  +  2  <  0 


2x2  +  2xy  +  y2  -  2x 

<  0 

m 

2  t  2 
x  +  xy  +  y  -  y 

1  A 

O 

These  inequalities  have  the  unique  solution  (22'"*  -  2^^  ,  2  -  2“^) 

3.  The  above  example  involves  three  inequalities.  It  seems  possible 
that  for  the  case  of  only  two  positive  definite  inequalities  rational 
solutions  will  always  exist.  The  answer  is  not  known  at  this  time. 
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